% 将 T_at_time 叠加228次，每次向上平移一行
time_step=700;
T_at_time= T(:, :, time_step);
T_at_time1 = T(:, :, time_step);  % 归一化温度场
clear new_image;
new_image = zeros(1001 + 228, 1001);

for i = 1:228
    new_image(i:i+size(T_at_time1, 1)-1, :) = new_image(i:i+size(T_at_time1, 1)-1, :) + T_at_time1;



end
new_image = new_image / max(new_image(:)) * 100;
% 显示最终结果

% 绘制最终的温度场
figure;
r_corrected = r * (2.27/0.0075)* (2.7/2);  % 横坐标 (r) 进行校正
z_corrected = z * (2.27/0.0075)* (2.7/2);;  % 纵坐标 (z) 进行校正

imagesc(r_corrected, z_corrected, new_image);  % 绘制二维温度场
% caxis([30 40]);

axis xy;  % 设置坐标轴方向正确
xlabel('r (mm)');
ylabel('z (mm)');
title(['Temperature Distribution at t = ', num2str(time_step*0.01), ' s']);
colorbar;  % 显示颜色条

% 设置 z 轴的范围
zlim([min(T_at_time(:)), max(T_at_time(:))]);  % 设置z轴的范围
% zlim([0, 100]);  % 归一化后最大温度为100
set(gca, 'ztick', []);  % 隐藏 z 轴刻度
set(gca, 'zticklabel', []);  % 隐藏 z 轴刻度标签

axis equal;  % 保证坐标轴比例相同

